(* Quadratura de f composta de Simpson,usando particao uniforme em n subintervalos de [a,b];
precisao casas decimais: prec *)

f[x_]:=Sin[x];
a=0;
b=Pi;
Nmax=100;
prec=10;

S[j_,m_]:=(b-a)/(6m)(f[(b-a) j/m ]+4 f[(b-a)(2j+1)/(2m) ]+f[(b-a)(j+1)/m ]);

Print["n valor quadratura"]; 
For[n=1,n≤Nmax,n++,
Print[n," ",
N[ Sum[S[i,n],{i,0,n-1}],prec]
]
];